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A COMPREHENSIVE  PROGRAM  FOR  AUTOMATING  RESIDUAL  STRESS  MEASUREMENTS 

BY  X-RAYS 

M.  R.  James  and  J.  B.  Cohen 
Northwestern  University 
Evanston,  Illinois  60201 

ABSTRACT 

The  software  design  of  a complete  computer  controlled  X-ray 
residual  stress  is  described.  The  program  will  operate  using  either 
a position  sensitive  detector  or  a diffractometer  for  data  accumulation. 

For  use  with  a diffractometer,  optimization  of  the  time  of  data 
accumulation  for  a prescribed  statistical  counting  error  is  presented. 

The  peak  location  routine, using  a multiple  pass  refinement  of  the  peak 
position,  is  detailed. 

The  software  allows  operator  intervention  during  operation  and 
specification  of  either  the  'two  tilt'  or  'sin3^'  technique.  The 
Marion-Cohen  technique  is  automatically  accomplished  if  oscillations 
in  d vs.  sin3*  are  evident.  Provisions  for  the  receiving  slit  movement, 
counting  statistical  analysis  and  instrumental  error  analysis  are 
provided.  The  operator  specifies  the  total  error  in  stress  he  will 

tolerate,  and  the  program  then  optimizes  the  counting  strategy. 
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CHAPTER  4 


4 . 1 INTRODUCTION 

Manual  measurement  of  residual  stress  requires  accumulating  X-ray 

counts  at  a minimum  of  three  2G  angles  over  the  peak  at  each  ^ inclination 

of  the  sample.  Long  collection  times  are  needed  to  accumulate  the  data 

making  the  procedure  tedious  and  time  consumming.  Hardware  controlled 

step  scanning  may  save  on  operator  time,  but  is  of  limited  value  since 

the  peak  position  and  breadth  varies  with  sample,  with  residual  stress 

and  with  psi  angle.  Indeed,  the  need  for  operator  control  has  been  the 

chief  reason  for  the  proliferation  of  the  'two  tilt'  method  over  the 
2 

'sin  f ' method,  although  more  precision  is  obtained  with  the  latter. 

Instrumentation  such  as  Fastress^ and  the  Shimadzu  X-Ray  Diffraction 

Stress  Analyzer^7^  offers  limited  automation  at  the  expense  of  precision 

(79) 

and  versatility.  Other  automated  systems'  7 require  the  collection  of 
many  data  points  in  the  peak  region  due  to  the  lack  of  on  line  control 
to  locate  the  peak.  None  of  these  instruments  offer  any  type  of  statis- 
tical analysis  or  optimization  of  data  accumulation. 

On  line  control  has  been  achieved  by  Kelly  and  Short and  Hayama 

/Q  1 \ 

and  Hashimoto'  by  using  step  scanning  techniques.  Both  use  counting 
statistics  to  estimate  the  statistical  error,  but  do  not  collect  data 
to  achieve  an  operator  specified  error  (i.e.,  they  do  not  optimize  the 
data  collection). 

Kelly  and  Short  designed  their  software  package  specifically  for 
the  ’two  tilt1  method  of  X-ray  stress  analysis  and  allow  only  a three 

*Ch.  1,  2 and  3 appeared  as  ONR  TR  No.  14  on  this  contract. 
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point  parabolic  fit  to  determine  the  peak  position.  Hayama  and  Hashimoto 
employed  the  parallel  beam  method  of  X-ray  optics  and  the  half-value 
breadth  method  to  define  the  diffraction  angle. 

The  system  to  be  described  includes  the  hardware  and  software  design 
for  computer  controlled  X-ray  measurement  of  residual  stress,  including 
sample  alignment,  optimized  data  accumulation  and  on  line  calculation 
of  the  stress  and  all  principal  errors.  The  system  is  designed  for  use 
with  either  a normal  diffractometer  with  step  scanning  or  with  a position 
sensitive  detector  (PSD)  to  collect  the  data.  The  operator  specifies 
the  error  in  stress  and  responds  to  certain  questions  on  the  equipment 
after  which  the  computer  devises  the  appropriate  counting  strategy. 

4 . 2 INSTRUMENTATION 

4.2.1  Hardware 

Control  of  the  diffractometer  was  executed  through  a PDP8/E  hardware/ 
software  system.  The  basic  DEC  PDP8/E  computer  having  16  K,  twelve-bit 
words  of  memory  and  dual  Dectape  transport  for  mass  storage  was  com- 
plimented with  interfacing  for  two  detectors,  a clock  to  serve  as  a time 

base,  and  a drive  for  three  SLO-SYN  stepping  motors.  The  basic  interface 

/ q9  \ 

arrangement  has  been  described  previous lyv  and  is  used  for  complete 
automation  of  the  Long-Term  X-Ray  Facility  at  Northwestern  University. 
rRecause  of  limited  funds,  sophisticated  items  such  as  absolute  angle 
decoders  or  addressable  axis  positions  and  counters  are  not  incorporated. 
Therefore,  all  motor  movements  and  counting  procedures  are  accomplished 
via  software  programing.  This  results  in  a versatile  system  that  can 
be  easily  adapted  (by  software  programming)  to  control  many  diverse 
experiments  in  the  field  of  X-ray  di f fraction. 1 
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4.2.2  Software  Package 

The  basis  of  the  software  system  is  a modified  version  of  OMSI-FOCAL 

tailored  specifically  to  meet  the  needs  of  the  X-ray  facility.  OMSI-FOCAL 

is  a Digital  Equipment  Corporation  FORTRAN-like  high  level  language 

(83 ) 

as  modified  by  the  State  of  Oregon  Museum  of  Science  and  Industry. 

User  defined  functions  have  been  created  bv  the  current  X-ray 
group  at  Northwestern  University  (including  the  present  author)  to 
allow; 

1)  Simultaneous  movement  of  three  stepping  motors. 

2)  Data  accumulation  by  fixed  time,  fixed  counts,  or  fixed 
monitor  counts. 

3)  Accumulation  of  data  from  two  detectors  simultaneously. 

4)  Two  hundred  ninety-six  directly  addressable  variables  and  992 
indirectly  addressable  arrayed  variables. 

5)  Accumulator  switching  to  enable  control  of  a teletype,  digital 
voltmeter,  thermae  controller,  and  a multichannel  analyzer  interface. 

This  function  allows  a change  in  the  operating  status  of  the  teletype 
(on-off)  or  a Thermae  controller  (from  a low  set  point  to  a high  set 
point),  monitoring  the  readout  of  a digital  voltmeter,  and  transfer  of 
data  from  a multichannel  analyzer  to  the  computer. 

This  version  of  Focal,  designated  by  the  X-ray  group  as  X-RAY 
FOCAL,  was  designed  to  alleviate  any  additional  machine  language  pro- 
gramming in  automating  X-ray  diffraction  experiments.  The  computer 
program,  STRESS,  to  automate  the  residual  stress  measurement  is  described 


in  Sec.  4.3. 
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4.3  PROGRAM  FEATURES 
4.3.1  Introduction 

The  following  features  are  implemented  in  the  automated  residual 
stress  package  called  STRESS. 

a.  Automatic  sample  alignment. 

2 

b.  Use  of  either  'two  tilt',  'sin  ^ ' or  ' Mar ion-Cohen 1 methods  of 
residual  stress  analysis  (see  Sec.  4.3.4). 

c.  Use  of  a step-scanning  diffractometer  or  a position  sensitive 
detector  system  for  data  accumulation. 

d.  Parafocusing  or  stationary  slit  'beam  optics'  (see  Sec.  4.3.5). 

e.  A least-squares  parabola  to  determine  the  peak  location  (see 
Sec.  4.3.6). 

f.  On-line  peak  location  with  operator  specified  accuracy  (see 
Sec.  4.3.7  and  4.3.8). 

g.  Specification  of  experimental  elastic  constants  or  single 
crystal  elastic  constants  (see  Sec.  4.3.9). 

h.  Instrumental  error  calculations  (see  Sec.  4.3.10). 

2 

i.  Plot  of  data  for  'sin  and  'Marion-Cohen'  methods. 

j.  Oscillation  of  sample  (see  Sec.  4.3.2). 

A master  flowchart  describing  the  overall  sequence  of  subprograms 
is  given  in  Fig.  4.1.  A thorough  description  of  each  feature  of  the 
program  follows.  The  master  dialog  sets  operating  parameters  for  the 
entire  residual  stress  package  and  is  explained  in  Sec.  4.3.2.  The 
program  may  branch  to  a sample  alignment  program  or  the  initial  dialogue 
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for  the  residual  stress  program  using  either  a step  scanning  diffractometer 
or  a position  sensitive  detector.  It  should  be  pointed  out  that,  as 
seen  in  Fig.  4.1,  if  the  stress  program  using  the  PSD,  STRPSD,  is  being 
used,  the  sample  alignment  program  cannot  be  implemented.  The  PSD 
needs  to  be  calibrated  to  change  linear  position  along  the  detector 
into  °20.  This  would  not  be  practical  in  determining  peak  positions  for 
multiple  peaks  because  an  absolute  20  along  the  PSD  would  not  be  known 
once  the  PSD  was  moved.  Thus,  the  PSD  is  only  practical  for  repetitive 
measurements  on  similar  samples  or  for  use  on  a dedicated  residual  stress 
device. 

Sample  input  and  output  are  provided  for  each  subprogram  with 
flowcharts  detailing  specific  sections.  Line  numbers  on  the  input-output 
figures  are  referred  to  in  the  text  for  convenience  to  the  reader. 


4.3.2  Master  Dialogue 

The  master  dialogue  sets  parameters  which  are  used  in  the  five 
subprograms  of  the  residual  stress  package  and  chains  to  the  appropriate 
subprogram  when  complete.  A sample  input  is  given  in  Fig.  4.2. 

The  master  program  chains  to  the  position  sensitive  detector  sub- 
program if  a positive  response  is  given  to  line  2 of  Fig.  4.2.  The 
sample  alignment  program  is  run  if  a positive  response  is  given  to  line 
9 in  Fig.  4.2. 

Two  options  affecting  the  stepping  motor  pulse  rate  are  necessary 
for  the  package  to  be  used  on  all  diffractometers  in  the  X-ray  facility 
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at  Northwestern.  By  responding  to  lines  3 and  4 in  Fig.  4.2  with  alphanumeric 
responses  of  a ' Y ' for  yes  or  a 'N'  for  no  the  program  8ets  the  proper 
pulse  speed  to  the  motors. 

Line  5 of  Fig.  4.2  determines  if  a solid  state  detector  is  to  be 
used.  Because  of  the  vibration  due  to  the  movement  of  liquid  , a 
delay  between  movement  of  the  detector  and  the  counting  of  photons  is 
made  to  prevent  the  detection  of  noise  generated  in  the  detector.  The 
delay  t ime  is  input  by  the  user  in  line  5,  Fig.  4.2.  The  deadtime  of 
the  counting  system  is  input  in  line  6. 

Specimen  oscillation  may  be  used  to  increase  the  number  of  crystals 
involved  in  diffraction  for  coarse  grained  materials.  This  method  is 
based  on  the  idea  that  the  stress  measured  by  X-ray  diffraction  is  the 
value  averaged  over  the  area  illuminated.  Small  oscillations  in  the  \jr 

2 

axis  will  not  greatly  affect  the  result  since  the  average  value  of  sin  i| / 
will  not  change  at  each  peak.  To  initiate  this  feature,  the  user  responds 
with  a * Y * to  line  7 of  Fig.  4.2  and  then  inputs  the  rocking  width  of 
the  oscillation.  The  axis  of  oscillation  is  the  axis. 

An  option  is  given  in  line  1 of  Fig.  4.2  to  skip  all  these  questions. 

This  default  scheme  assumes  the  use  of  the  Picker  diffractometer  with  a 
scintillation  detector,  and  sets  the  deadtime  to  1.5  u,sec.  The  program 
chains  directly  to  the  initial  dialogue  for  the  residual  stress  program  using 
a diffractometer  or  to  the  sample  alignment  program  if  requested. 

4.3.3  Sample  Alignment 

The  displacement  of  the  sample  from  the  true  center  of  the  diffractometer 
can  be  estimated  from  the  slope  of  the  lattice  parameter  vs.  Nelson-Riley 
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factor,  A subprogram,  ALIGN,  automatically  locates 

the  Bragg  angle  of  any  peak  by  defining  the  peak  position  from  the  apex 
of  a parabola  fit  to  the  top  15  pet  of  the  diffraction  profile.  After 
determining  three  or  more  peaks  specified  by  the  operator,  the  lattice 
parameter  and  Nelson-Riley  factor  are  determined.  The  extrapolated 
lattice  parameter  is  determined  from  a linear  least-squares  fit  to  the 
lattice  parameter  vs.  Nelson-Riley  factor.  The  sample  displacement  is 
determined  from  the  slope  of  the  least  squares  fit  after  which  the 
operator  may  manually  or  automatically  reposition  the  sample  and  repeat 
the  alignment  or  chain  to  the  initial  dialogue  for  the  diffractometer 
version  of  the  automated  residual  stress  program. 

Figure  4.3  reproduces  the  input  and  output  of  the  alignment 
program.  Only  three  data  points  are  used  in  the  parabolic  fit.  The 
operator  needs  only  specify  the  current  angular  settings,  the  initial 
20  position  from  where  to  begin  the  first  step  scan  to  find  the  rough 
peak  location  and  the  20  step  increments  for  the  first  and  second  step 
scans  (designated  DELTA  1 and  DELTA  2,  respectively)  to  refine  the  peak 
position.  This  peak  location  scheme  is  identical  to  that  used  in  the 
residual  stress  program,  a detailed  description  of  the  procedure  being 
given  in  Sec.  4.3.7.  The  only  change  is  that  the  operator  has  the  option 
in  the  alignment  program  to  specify  the  accuracy  of  each  peak  location 
(in  °20).  If  he  does  not,  a fixed  count  of  10,000  is  obtained 

for  each  data  point. 

The  peak  location  and  the  error  on  the  peak  location  (one  standard 
deviation)  are  output.  For  convenience  the  peak  intensity  and  the  step 
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increment  between  the  final  data  points,  DELTA,  are  printed.  If  the 
operator  desires  to  repeat  the  alignment;  e.g.,  if  he  repositions  the 
sample  and  wants  to  check  it,  this  information  enables  him  to  judge 
more  accurately  the  precision  of  peak  location  and  the  two  step  increments. 

The  extrapolated  lattice  parameter,  slope  and  correlation  coefficient 
from  the  lattice  parameter  vs.  Nelson-Riley  plot  are  printed.  The  sample 
displacement  is  typed  after  the  operator  gives 

the  goniometer  radius.  The  alignment  program  may  be  repeated  or  the 
initial  dialog  program  for  the  diffractometer  version  for  residual 
stress  determination  may  be  run. 

4.3.4  Method  of  Residual  Stress  Analysis 

2 

The  'two  tilt',  'sin  ty'  or  'Marion-Cohen ' methods  of  residual  stress 

analysis  are  available  by  responding  with  a * Y ' for  'yes'  to  the  appro- 

★ 

priate  question  (lines  1-3  in  Fig.  4.4  ).  The  specifics  of  each  method 
are  detailed  below.  This  applies  to  both  the  diffractometer  and  position 
sensitive  detector  version. 

(A)  'Two  Tilt’  Method 

If  the  operator  chooses  this  method  only  two  tilts  of  the  specimen 
are  employed  to  determine  the  residual  stress.  The  * angles  are 
requested  for  each  tilt  (line  19-20,  Fig.  4.5).  The  stress  is  determined 
using  either  an  exact  formalization  or  an  approximate 

equation  relating  the  stress  directly  to  the  peak  shift  through  a stress 
constant,  K.  The  choice  is  made  when  the  operator  inputs  the  elastic 

k 

Because  of  the  many  options,  two  input  dialogs  are  given  in  Fig.  4.4  and 
Fig.  4.5.  Line  numbers  and  the  figure  number  are  given  wherever  a refer- 
ence to  each  figure  is  given.  Typical  output  is  given  in  Fig.  4.6  for  the 
two  tilt  and  Fig.  4.7  for  the  sin2*  and  Marion-Cohen  method. 
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FIGURE  A. A Sample  input  for  diffractometer  version  of  STRESS. 

Data  is  input  for  a typical  measurement  using  the 
sin^f  technique.  Copyright  Northwestern  University,  1976. 
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FIGURE  4.5  Sample  Input  for  diffractometer  version  of 

STRESS.  Data  Input  for  a typical  measurement 
using  the  two-tilt  technique. 

Copyright  Northwestern  University,  1976. 
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S'PESS  121.27  MPA  < 14653  °Sl  ) 

STA' I ST  I CAL  EDD0P  ♦/-  9.06  MPA  ( ♦/-  1319  5?I  ) 


FIGURE  4.6  Sample  output  for  STRESS  package  using  the 
two-tilt  technique.  Background  feature  was 
used  and  scattering  facter  correction  made  to 
the  raw  data. 

Copyright  Northwestern  University,  1976. 
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FIGURE  4.7  Sample  output  for  STRESS  package  using  the  sin  f method 
with  Marlon-Cohen  texture  analysis. 
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constants  (lines  21-23,  Fig.  4.5). 


The  statistical  counting  error  in  the 
stress  for  the  two  tilt  method  is  calculated 

depending  on  whether  an  exact  or  approximate  formalization  is 

used. 


(B)  Sin  \li  Method 

2 

Up  to  eight  tilts  are  allowed  when  using  the  sin  ^ method,  as 
shown  in  lines  17-20,  Fig.  4.4.  The  stress  is  calculated 


and  the  statistical  counting  error  in  stress  determined. 

/ Q C \ 

A linear  least  squares  method  given  by  Ezekiel  and  Fox  is  utilized 

2 

to  calculate  the  slope  of  d vs.  sin  The  intercept,  correlation 

if 

coefficient  and  standard  deviation  of  the  slope  are  calculated  and 

2 

printed.  A plot  of  d^  vs.  sin  ^ is  given  if  requested. 


(C)  Marion-Cohen  Method 

The  use  of  the  Marion-Cohen  method  to  enable  the  determination  of 
the  macroscopic  residual  stress  in  the  presence  of  oscillations  in  d 

* 

2 2 
vs.  sin  ijr  has  been  presented.  All  features  of  the  sin  f 

2 

method  are  incorporated  in  this  technique.  After  the  normal  sin 

analysis,  the  correlation  coef f icientv  is  tested  to  judge  if  the  fit 
2 

of  d vs.  sin  <ii  is  linear  within  a 90  pet  confidence  level  or  if  oscil- 
* 

lations  may  be  present.  If  the  latter  is  true,  the  Marion-Cohen  analysis 
is  implemented. 

The  required  distribution  or  texture  function  is 

determined  using  the  peak  intensity  at  each  ^ tilt.  Both  the  maximum 
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and  background  Intensity  are  obtained  with  the  receiving  slit  positioned 

on  the  goniometer  circle.  This  is  to  insure  that  the  same  size  arc  on 

the  Debye  cone  is  recorded  for  each  tilt  which  allows  for  normalization 

(78) 

of  the  distribution  function.  The  background  is  subtracted  from  the 

•k 

maximum  intensity  and  the  net  value  corrected  for  absorption  as  described 
by  Marion.  These  values  are  normalized  with  the  maximum  value  as 

unity.  A multiple  least-squares  fit  is  performed  and 

(d  -d  ),  d , and  the  stress  are  calculated.  A typical  output  is  given  in  Fig.  4.7. 
max  fi  B 

4.3.5  Focusing  Technique 

The  operator  is  queried  whether  or  not  he  wishes  to  have  the 
receiving  slit  moved  in  line  4,  Fig.  4.4.  This  question  effectively 
determines  if  parafocusing  geometry  (movement  of  the  receiving  slit) 
or  non-focusing  geometry  is  implemented.  An  approximate  peak  position 
given  by  the  operator  in  the  initial  dialog,  line  6 of  Fig.  4.5,  is 
used  in  the  calculation  of  the  distance  of  movement. 

Instrumental  errors  are  calculated  using  the  appropriate  sample  to 
receiving  slit  distance  determined  for  each  technique. 

4.3.6  Number  of  Data  Points  in  Peak  Location 

Independent  of  the  method  of  X-ray  stress  analysis,  the  operator 
specifies  the  number  of  observation  positions  for  the  parabolic  fit  in 
line  14  of  Fig.  4.4.  At  least  three  data  points  must  be  used  to  fit  a 
parabola  to  the  profile.  A least-squares  routine  is 

used  to  locate  the  apex  of  the  parabola  and  requires  only  an  odd  number 
of  data  points. 

★ 

The  absorption  correction  is  necessary  because  of  the  change  in  the  area 
of  the  specimen  irradiated  in  going  from  one  ^ angle  to  another  and  the 
relative  distances  travelled  by  the  Incident  and  diffracted  rays  through 
the  specimen. 
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4.3.7  On-Line  Peak  Location 

Automated  location  of  the  peak  of  a diffraction  profile  requires 
first  determining  the  region  of  fit  of  a smooth  curve  to  the  profile, 
obtaining  the  data  at  the  proper  angular  positions  to  the  desired  pre- 
cision and  processing  the  data.  A flowchart  covering  the  peak  location 
is  given  in  Fig.  4.8. 

The  region  of  fit  for  a parabola  is  chosen  using  the  'top  15  pet 
rule'.  An  option  is  supplied  (line  10,  Fig.  4.5)  allowing  the  back- 

ground to  be  subtracted  prior  to  determining  the  region. 

The  angular  position  at  which 

the  background  intensity  is  to  be  obtained  must  be  input.  A constant 
value  for  the  background  is  used  because  the  correction  is  only  made  in 
determining  the  region  of  parabolic  fit  and  is  not  crucial. 

An  approximate  position  for  the  maximum  in  intensity  is  found 
utilizing  a step  scan  initiated  from  an  angular  position  (input  in  line 
6,  Fig.  4.5)  below  the  peak.  A large  step  increment,  specified  by  the 
operator  in  line  7 of  Fig.  4.5,  may  be  used  to  speed  up  this  first  scan. 
A fixed  count  of  1000  is  used  unless  the  operator  desires  to  decrease 
the  preset  counts  if  a weak  peak  is  being  examined  (line  11,  Fig.  4.5). 
The  maximum  intensity  is  stored  with  the  scan  continuing  until  a 
count  rate  of  less  than  90  pet  of  the  maximum  is  found.  This  allows  a 
99  pet  (three  sigma)  confidence  in  the  statistical  count  rate  (if  1000 
counts  are  used)  to  be  sure  the  maximum  has  been  passed.  With  the 
maximum  count  rate  and  its  angular  position  stored,  two  step  scans  down 
each  side  of  the  peak  in  smaller  step  increments  are  initiated.  The 
increment  for  this  second  step  scan  to  further  refine  the  peak  is 
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FIGURE  4-8  Flowchart  for  peak  location 
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specified  in  line  9,  Fig.  4.5.  The  maximum  count  rate  is  adjusted  if 
a greater  intensity  than  heretofore  assigned  to  the  maximum  is  obtained 
at  any  position  during  the  second  step  scan  using  the  smaller  step 
increments.  The  step  scans  continue  until  angular  positions,  20^  and 
20^,  corresponding  to  85  pet  of  the  maximum  intensity  (with  background 
subtracted  if  desired)  are  located.  These  two  positions,  20^  and  20^, 
fix  the  beginning  and  end  of  the  region  of  the  parabolic  fit. 

(A)  Data  Accumulation 

Using  the  two  end  points,  20^  and  2 0^,  and  the  working  origin  as 
half  way  between,  a preliminary  three  point  parabolic  fit,  using  fixed 
counts  of  5000,  is  made.  All  angular  intervals  are  rounded  to  the 
nearest  .005°20  when  they  are  calculated  because  the  20  SLO-SYN  motor 
can  only  be  positioned  to  this  accuracy.  The  data  is  corrected  for 
Lorentz-polarization,  absorption  and  scattering  factors 
and  a parabola  fit  to  the  data  is  obtained.  Based  on 

this  parabola  the  angular  positions  having  85  pet  of  the  peak  intensity 
are  determined.  The  step  increment  is  determined  using  the  desired 
number  of  data  points  (see  Sec.  4.3.6).  The  scan  is  repeated  once  more 
using  the  desired  number  of  observation  points  to  obtain  data  on  which 
to  estimate  the  necessary  preset  counts  to  achieve  the  desired  accuracy 
(see  Sec.  4.3.8).  The  final  data  scan  i6  made  using  the  same  angular 
positions  as  the  previous  scan  but  for  the  calculated  preset  counts. 

Based  on  numerous  preliminary  tests,  this  method  of  five  refining 
step  scans  was  judged  to  be  the  fastest  and  most  reproducible  method  of 
determining  the  peak  position,  especially  for  broad  peaks.  In  summary, 
the  five  scans  are: 
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1)  Find  approximate  peak  using  large  step  increments.  This  normally 
takes  about  30  seconds,  depending  on  the  intensity. 

2)  Find  region  of  fit  using  small  step  increments.  For  a sharp 
peak  this  takes  approximately  30  seconds  and  60  seconds  for  a broad  peak. 

3)  Refine  peak  location  using  parabolic  fit  to  three  angular 
positions  in  region  of  fit  (»».  60  seconds). 

4)  Obtain  data  rapidly  at  final  angular  positions  to  determine 
necessary  preset  counts  (~  30  seconds), 

5)  Obtain  final  data. 

(B)  Processing  of  Data 

The  raw  data  is  corrected  for  the  intensity  correction  factors 

for  the  Lorentz-polarization  absorption  factors,  respectively. 
Correction  for  the  atomic  scattering  factor  is  impieme.ited  when  desired 
by  the  operator.  The  operator  must  specify  two  values  of  the  scattering 
factor  at  their  respective  value  of  sinO/J  as  indicated  in  lines  17-18, 

Fig.  4.5.  A straight  line  is  fit  to  this  data  and  the  slope  used  to 
calculate  the  appropriate  scattering  factor  at  each  angular  position. 

The  data  is  d ivided  by  the  Lorentz-polarization,  absorption  and  scat- 
tering factor  squared.  The  temperature  factor  is  not  included  because 
it  is  so  insensitive  to  20  but  this  would  be  included 

if  the  operator  inputs  fe  when  asked  for  f. 

Peak  location  is  determined  with  the  counting 

statistical  error  in  the  peak  calculated  for  fixed  count 

data  accumulation.  A report  of  the  peak  angular  position,  intensity  at 
each  observation  point,  the  appropriate  correction  factor,  the  time  of  data 
accumulation  at  each  point  and  the  statistical  counting  error  is  typed. 
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This  procedure  of  peak  location  and  output  is  continued  until  all 

ilf  tilts  are  complete.  Finally,  the  residual  stress,  counting  error  and 
★ 

instrumental  error  are  calculated  and  a report  typed.  A typical  output 

2 

is  given  in  Fig.  4.6  for  the  two  tilt  method  and  Fig.  4.7  for  the  sin  -fr 
and  Marion-Cohen  methods. 

4.3.8  Optimization  of  Data  Collection  for  a Specified  Accuracy 

In  order  to  determine  the  appropriate  counting  time,  the  program 
must  be  told  the  desired  accuracy  in  the  peak  position.  The  program 
allows  the  error  to  be  input  in  the  initial  dialog  (line  29,  Fig.  4.4) 
in  terms  of  °20  for  the  peak  location  or  in  terms  of  stress  from  which 
the  appropriate  error  in  the  location  of  each  peak  is  determined. 

In  the  first  case,  direct  input  of  the  accuracy  for  the  peak  location, 
the  operator  specifies  in  °20  the  precision  to  which  each  peak  should  be 
located.  One  value  is  input  and  each  peak  located  to  that  accuracy  using 
only  the  counting  statistics.  This  is  shown  in  lines  24-25,  Fig.  4.5. 

If  the  total  error  is  given  in  terms  of  stress,  as  in  line  29,  Fig. 

4.4,  the  computer  must  back  calculate  the  error  in  peak  location  as  follows 

a)  Estimate  the  instrumental  error.  This  is  done  using  the 
approximate  peak  location  as  specified  by  the  operator  and  other  data 
given  by  the  operator  (see  Sec . 4 . 3 . 10  for  discussion  of  instrumental 
error  input). 

b)  The  error,  E,  in  stress  due  to  counting  statistics  alone  is 

then  calculated  from:  , 

r 2 2 1 

E(counting  .t.tl.'ic.)  " LE  <T0T4L'-E  (INSTRUMENTAL  ERROR)  | 

(4.1) 

The  exact  instrumental  error  and  the  correct  sign  is  printed.  This 
may  be  subtracted  from  the  measured  stress  (if  the  input  parameters  are 
correct)  or  may  treated  simply  as  an  additional  error.  Eq.  4.1  may  be 
used  to  obtain  a total  error  if  the  latter  treatment  is  preferred. 
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c)  The  error  in  peak  location  is  then  back  calculated  from  the 
counting  statistical  error  (which  ia  in  units  of  stress)  using  equations 
pertaining  to  the  choice  of  the  measurement. 

The  calculations 

are  made  to  give  the  same  statistical  error  at  each  \ji  tilt.  The  computer 
types  out  the  error  in  the  peak  location  to  be  sure  it  is  acceptable 
to  the  operator.  This  error  represents  one  standard  deviation  confidence 
limits  and  all  peaks  are  determined  to  this  precision. 

The  actual  determination  of  the  preset  ccwnts  takes  place  in  the 
fourth  step  scan  (see  Sec.  4.3.7).  From  the  third  step  scan,  the  final 
angular  positions  are  determined.  In  the  fourth  step  scan,  data  is 
taken  rapidly  at  these  angular  positions.  Using  this  data  the  sumnations 
necessary  to  derive  the  error  in  peak  location  for  fixed 

counts  can  be  calculated.  Setting  the  specified  counting  error  equal  to 

2 

n (2  tip)  enables  the  minimum  preset  counts  to  obtain  the 

desired  accuracy  to  be  calculated.  The  final  data  is  taken  for  the 
necessary  fixed  counts. 

This  process  optimizes  the  time  of  data  collection  for  a specified 
accuracy.  The  calculation  time  is  minimal  and,  because  all  but  the  fourth 
step  scan  are  necessary  to  determine  the  position  of  the  observation 
points,  very  little  extra  time  is  taken  to  optimize  the  final  data 
accumulation.  The  fourth  step  scan,  usually  taking  approximately  30 
sec.,  is  necessary  to  obtain  an  accurate  prediction  of  the  summation  term. 

Preliminary  tests  showed  that  if  the  data  obtained  in  the 
third  step  scan  is  used  to  back  calculate  the  fixed  counts 
the  estimate  of  the  preset  counts  is  not  accurate  because  the  angular 
positions  of  the  data  points  are  different  between  the  third  and  final 


step  scans. 
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If  the  operator  requests  an  estimate  of  the  total  time  (line  18, 

Fig.  A. 5)  the  program  will  print  the  total  length  of  time  of  the  measurement 
based  on  the  time  necessary  to  obtain  the  first  peak.  At  this  point,  the 
operator  may  change  the  counting  error  if  the  time  is  not  acceptable 
from  which  a new  total  time  of  analysis,  based  on  the  current  data,  is 
calculated  by  the  computer  and  printed  out.  This  sequence  is  repeated 
until  the  operator  accepts  the  total  time  and  operation  continues. 

4.3.9  Specification  of  the  X-Ray  Elastic  Constants 

The  operator  may  specify  the  elastic  constants  in.  one  of  three 
ways  and  may  use  units  of  either  MPa  or  PS1  as  shown  in  line  21  of 
Fig.  4.5. 

If  the  approximate  formula  for  the  two  tilt  method  is  to 

be  used  the  stress  constant,  K,  must  be  input.  The  stress  constant 
depends  on  the  ^ inclination  and  it  is  up  to  the  operator  to  use  the 
proper  tilt. 

The  best  method  to  input  the  elastic  constants  is  to  use  measured 
values  from  the  literature.  These  must  be  input  as  (l+vl/Eand  may  be 
determined  from  mechanically  measured  bulk  elastic  constants  or  prefer- 
rably  from  values  measured  by  X-ray  diffraction.  A good  list  is  given 
in  ref.  90. 

The  single  crystal  elastic  constants  may  be  used  to  calculate  the 
elastic  constants  of  a polycrystalline  material  (see  lines  24-27, 

Fig.  4.4),  The  theoretical  methods  vary  according  to  the  fundamental 
boundary  conditions  used  to  account  for  the  influence  of  elastic  aniso- 
trophy  and  grain  interactions.  A satisfying  solution  applicable  to  the 
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/ q7  \ 

X-ray  elastic  constants  has  not  yet  been  achieved.  Modified 

approximations  of  Voight^7^  and  Ruess^75^  for  the  calculation  of  elastic 
constants  for  quasi-isotropic  polycrystals  are  employed.  The  equations 
for  the  Voight  model  are: 


.VOIGHT 

-v/£  - 


S11<2SU  + 2S12 


S44>  + S12(3S44 


4S12> 


2S44  + 6<S11 


S12> 


(4.2) 


and 


1+v  _ _ VOIGHT  _ 5 (S11  " S12)  S44 

E ^ 2}X-RAY  2S44  + 6(SU  - S12) 


where  are  the  elastic-compliance  constants  for  the  corresponding 
single  crystals.  The  analysis  by  Ruess  yields; 
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As  pointed  out  by  Neerfield^^  and  Hill^^  values  of  and  S^/2 
usually  lie  between  the  values  calculated  by  the  limiting  assumptions 
of  Voight  and  Ruess.  Therefore,  arithmetic  mean  values  of  the  elastic 
constants  calculated  from  the  two  models  are  used  in  the  computer  progran . 
Input  of  the  three  single  crystal  constants,  S^.  an<^  are 

requested  in  lines  24-26  of  Fig.  4.4.  The  calculated  value  of  1+v/E 
is  typed  and  used  to  determine  the  stress. 

4.3.10  Calculation  of  Instrumental  Error 

The  instrumental  error  is  determined  if  requested.  The  necessary 
information  input  in  lines  30-33  of  Fig.  4.4  includes  the  divergent 
slit,  the  sample  displacement,  the  ty-axis  missetting  and 

the  goniometer  radius,  if  these  parameters  are  not  input,  the  instrumental 
error  is  not  printed  in  the  output,  as  in  Fig.  4.6.  The  exact 
aberration  (and  its  sign)  due  to  horlzantal  beam  divergence, 
sample  displacement,  and  ^-axis  displacement, 

are  determined  after  all  peak  positions  have  been  determined. 

In  the  initial  dialog,  the  operator  may  desire  a preset  total  error 
in  stress.  From  this,  the  statistical  counting  error  must  be  obtained 
using  Eq.  4.1.  Therefore,  an  estimate  of  the  instrumental  error  is 
made  by  using  the  approximate  peak  location  input  by  the  operator  at 
line  6,  Fig.  4.5.  This  gives  a reasonable  estimate  of  the  instrumental 
error  before  the  actual  peak  angles  are  determined  in  order  to  calculate 
the  desired  counting  error. 

4.3.11  Position  Sensitive  Detector  Program 

The  use  of  a position  sensitive  detector  (PSD1  for  residual  stress 
measurements  has  been  treated  previously.  A complete  description  of 
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the  detector  system  and  performance  was  presented  In  ONR  TR  Nos.  11,  12 
and  15.  The  PSD  eliminates  all  2b  motion  and  enables  data  over  the 
entire  diffraction  profile  to  be  accumulated  simultaneously.  The  data 
is  stored  in  a multichannel  analyzer  (MCA). 

The  automated  residual  stress  package  allows  the  data  to  be 
transferred  from  the  MCA  to  the  computer  for  processing.  The  program 
features  are  described  below  and  assume  the  reader  has  sufficient 
knowledge  of  operation  of  the  PSD  system.  Also,  previous  calibration 
of  the  detector  is  necessary. 

Figure  4.9  gives  a typical  dialogue  for  the  position  sensitive 

2 

detector  program.  The  two  tilt,  sin  ^ or  Marion-Cohen  methods  may  be 
used  with  the  same  features  as  explained  in  Sec.  4.3.2.  In  order  to 
transfer  the  data  from  the  MCA  to  the  computer,  the  first  channel  and 
last  channel  on  the  MCA  to  be  transferred  must  be  input  (lines  2-3, 

Fig.  4.9). 

Information  to  determine  the  calibration  of  position  along  the 
detector  to  °20  must  be  input  in  lines  4-6  of  Fig.  4.9.  This  data 
must  be  determined  from  a calibration  sample. 

The  linear  calibration  constant  in  °20/channel  on  the  MCA  is  input  in 
line  4.  The  absolute  reference  angle  in  °20  is  input  with  the 
corresponding  channel  on  the  MCA. 

The  Lorentz-polarization,  absorption  and  scattering  factor  corrections 
are  made  to  the  data  after  the  angular  position  of  each  data  point  is 
determined  by  the  computer  program.  The  program  next  determines  the 
region  of  parabolic  fit  by  including  all  data  points  within  85  pet  of 
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FIGURE  4.9  Input/output  from  PSD  version  of  STRESS 
Copyright  Northwestern  University,  1976. 
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the  maximum  intensity.  A least-squares  parabolic  fit  is 

then  made  utilizing  all  the  data  points  within  the  region.  Statistical 

counting  errors  are  determined  for  the  peak  location. 

The  number  of  ^ tilts,  elastic  constant  data  and  instrumental  error 
data  are  similar  in  form  to  the  program  section  for  the  diffractometer. 
All  calculations  for  the  stress,  counting  errors  and  instrumental  errors 
are  identical  to  those  previously  described  in  the  diffractometer  version 
of  the  residual  stress  package. 

Due  to  the  physical  size  of  the  PSD,  it  is  not  possible  to  place 
the  detector  on  the  focusing  circle  in  back  reflection  geometry. 

This  is  clearly  seen  in  Fig.  4.10  where  the  perpendicular  bisector  of 
the  detector  anode  is  placed  at  156°20,  a typical  back  reflection 
angle  for  use  with  steel  specimens.  The  length  of  the  anode  and  the 
surrounding  packaging  require  the  detector  to  be  placed  behind  the  X-ray 
tube  and  off  the  focusing  circle  by  5 cm. 

For  the  portable  residual  stress  analyzer  to  be  described  in  ONR 
TR  No.  18  this  defocussing  is  minimized  because  of  the  size  of  the  X-ray 


tube. 
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J 

FIGURE  4.10  Top  view  of  Picker  diffractometer  showing  the  location  of 
the  position  sensitive  detector  (A'.  Due  to  the  physical 
size  of  the  detector,  it  was  placed  3 cm  from  the  focusing 
circle  at  ty=(P  . 
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